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ABSTRACT 

The aromatic features in M101 were studied spectroscopically and photometrically using observa- 
tions from all three instruments on the Spitzer Space Telescope. The global SED of M101 shows 
strong aromatic feature (commonly called PAH feature) emission. The spatially resolved spectral and 
photometric measurements of the aromatic feature emission show strong variations with significantly 
weaker emission at larger radii. We compare these variations with changes in the ionization index 
(as measured by [Nelll] / [Nell] and [SIV/SIII], which we probe over the ranges 0.03-20 and 0.044-15 
respectively) and metallicity (expressed as log(0/H) + 12, which ranges from 8.1 to 8.8). Over these 
ranges, the spectroscopic equivalent widths of the aromatic features from 7 HII regions and the nu- 
cleus were found to correlate better with ionization index than with metallicity. This implies that 
the weakening of the aromatic emission in massive star forming regions is due primarily to processing 
of the dust grains in these environments, not to differences in how they form (although formation 
could still be important on a secondary basis). The behavior of the correlation between the aromatic 
feature equivalent widths and ionization index can be described as a constant equivalent width until 
a threshold in ionization index is reached ([NeIII]/[NcII] ~ 1), above which the equivalent widths 
decrease with a power law dependence. This behavior for M101 HII regions is a l so see n for the sample 
of starburst galaxies presented in the companion study of lEngelbracht et al.l (|2008h which expands 
the range of [NeIII]/[NeII] ratios to 0.03-25 and log(0/H)+12 values to 7.1-8.8. The form of the 
correlation explains seemingly contradictory results present in the literature. The behavior of the 
ratios of different aromatic features versus ionization index does not follow the predictions of existing 
PAH models of the aromatic features implying a more complex origin of the aromatic emission in 
massive star forming regions. 

Subject headings: galaxies: individual (M101) — galaxies: spiral — galaxies: ISM — dust, extinction 



1. INTRODUCTION 

The interpretation of infrared measurements of star 
forming galaxies rests on an understanding of the proper- 
ties of the dust producing the emission. In addition, the 
dependence of dust properties on environment provides 
important clues to the dust composition and formation. 
An important constituent of the interstellar dust is large, 
aromatic molecules/grains whose emission often domi- 
nat es the mid-infrared output of HII regions and galax- 
ies (IRoche et al.l 119911 : IMadden et ail 120061 : ISmith et all 
I2007bh . 

The aromatic features were first discovered in observa- 
tions of planetary nebulae where a feature at 11.3 /im 
was broader than the observed atomic emission lines 
(|Gillett et al.| [T973). The number of these features has 
grown to include man y features wi t h wa v elengths be- 
tween 3 and 18 u m (jWerner et all 12004 ISmith et all 
120041: lTielensl[2005T) . The strongest aromatics are seen 
at 3.3, 6.2, 7.7, 8.6, 11.3, 12.7, and 17.1 fim and all 
the aromatics have been identified with C-H and C- 
C bending and stretchi ng modes of h ydrocarbons con- 
taining aromatic ring s (|Tielensl [20051) . The shape and 
strength of the aromatics have been seen to vary in sin- 
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gle sources (IWerner et all 12 004V a mong v arious Galac- 
tic sourc es (Ivan Diedenhoven et all [2004D. and anion; 
galaxies (lEngelbracht et al.l 120051: IMadden et all 1200 
IWu et aljl2003 r~ 

A number of different materials have been proposed 
as the carriers of these features, including Hydro- 
genat ed Amorphous Carbon (HAC, iDulev fe Wi lliams 
1983ft . Quenched Carbonaceous Composites (QCC 



Sakata et al.lll984T), Polvcyclic A romati c Hydrocarbons 
(PAH s, lAllamandola et al 119851). coal dPapoular et al.l 



I1989T ) . and nanodiamonds ( Jones fe d'HendecourtJl200*fj[ r 
The leading candidate material for the aromatics is PAH 
molecules and they are often included in the modi- 
fied "astronomical PAHs" form in dust grain models 
(|Li fe Draindl200lHZubko et alJl2004f ). We refer to these 
features as the aromatic features to allow for discussions 
of the observations without making an assumption about 
the material that is responsible fo r them . We refer the 
interested reader to iPeeters et al.l ([2004) for a detailed 
discussion of these features. 

The observed strengths of the aromatic features do 
not vary much in most normal l uminosity galaxies 
with app r oxima t ely solar met a llicities ([Roche et al J 19911 
ILu et al l 1200! ISmith et~all l2007br ) or massive star- 
bursts ([Brandl et al.l I2006T) . The strength of these 
featu res do weaken significantly in active galactic nu- 
clei (IRoche et al.l 119911: iLutz et al.l 119981: ISmith et all 
l2007bh and star bursting galaxies with fairly low 
metallicities ([Roche et al.lll991b lEngelbracht et alJl2005t 
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Galliano et alll2005t iMadden et al .1120061 : 1 Wu et alJl20M 
Engelbracht et aiT ~ l2008h . The variations seen in ac- 
tive galactic nuclei are not considered for this paper, 
which only discusses those clearly associated with nearby 
massive star formation (e.g., starburst galaxies and HII 
regions). Observations with the Infrared Space Ob- 
servatory (ISO) found that the strength of the aro- 
matic features was correlated with a galaxy's metallic- 
ity (probed by log(0/H)+12) and ionization (probed by 
[Nelll] / [Nell] ) , but were unable to determ ine which was 
the dominant effect (jMadden et al J 12006). The correla- 
tions are such that the aromatic features are weakest for 
low metallicities and high ionizations. 

If the correlation is better with metallicity than ioniza- 
tion, then this implies that the formation of the aromat- 
ics (from the constituent atoms in circumstellar or inter- 
stellar environments) happens less efficiently than that of 
small grains as the metallicity decreases. On the other 
hand, if the better correlation is with ionization, then this 
implies that there is increased processing of the aromatic 
feature carriers (modification and/or destruction in the 
immediate environment) as the radiation field becomes 
harder. Throughout this paper, we use the metallicity 
as a probe of the formation of aromatics from the con- 
stitute atoms and ionization as a probe of the process- 
ing of the aromatics in the immediate environment. The 
metallicity and ionization in massive star forming regions 
is roughly correlated making it nontrivial to determine 
which is the dominate correlation. 

Observations with Spitzer have refined the observa- 
tional picture of the variations of the aromatic features 
among galaxies, but have yet to provide a clear pic- 
ture of the origin of the variations. Evidence that the 
aromatic features are absent at the lowest metallicities 
probed by starburst galaxies was provided by the IRS 
(Infrared Spectrograph) spectrum of SBS 033 5-052 which 
has alog(0/H)+12=7.2 (jHouck et al.ll2004af ) . The weak- 
ness/absence of aromatic features at metallicities below 
log(Q/H) +12=8.2 in a sample of starburst galaxies was 
shown by Engelbra cht et al"T(|2005h using IRAC (Infrared 
Camera) and MIPS (Multiband Imaging Photometry for 
Spitzer) phot ometry confirme d in a few cases with IRS 
spectroscopy. IWu et all (|2006l ) studied a small sample of 
Blue Compact Dwarf galaxies (including SBS 0335-052) 
with IRS spectroscopy and found similar correlations be- 
tween aromatic feature equiv alent width a n d met allicity 
or radiation field hardness. iBrandl et al.l (|2006f ) found 
that the equivalent width of the aromatic features did 
not vary significantly over a factor of 10 in [NeIII]/[NeII] 
for a sample of starburst galaxies. Using the HII classi- 
fied galaxies in the SINGS (Spitz er Nearby Galaxies Sur- 
vey) sample [Smith et ahl (|2007bf ) found variations in the 
summed strength of the aromatic features when ratioed 
to the total IR flux with some indication of a threshold at 
log(0/H)+12=8.1. The metallicity thresho ld in the aro- 
matic feature strengths was confirmed by iDraine et al.l 
(|2007h where the broadband spectra energy distributions 
(SEDs) of the SINGS galaxies were modeled. These 
models indicated a fairly abrupt reduction by a factor 
of three occurs in the fraction of PAHs near a metallicity 
of log(0/H)+12 = 8.1. 

All this variation among galaxies raises the question: 
do the aromatics vary significantly among HII regions 
in a single galaxy? The study of the simpler environ- 



ments of HII regions may also present a clearer pic- 
ture of the cause of the aromatic feature variations 
as HII regions have much more homogeneous physical 
conditions than whole galaxies. The ideal galaxy to 
answer this question is the giant spiral galaxy M101. 
M101 is a large (> 30'diameter), face- on spiral galaxy 
(SAB cd) at a distance of 6.7 Mpc (|Freedman et al.l 
l200l . It has one of the largest metallicity gradients 
known with metallicities [log(0/H) + 12] ranging from 
8.8 in the nuc leus to 7.4 in an HII region 41 kpc 
from the center (IZaritskv et alJll994Hvan Zee et alj[l99a 
iKennicutt et all l2003al IBresolinl l2006h . In addition, 
there are high quality electron temperature and metal- 
licity determinations for a samp le of M101 HII regions 
with metallicit ies from 7.4 to 8.8 (|Kennicutt et al.l F2003a: 
IBresolinl 12006). This large apparent size, the numerous 
HII regions, the large metallicity gradient, and the high 
quality metallicity measurements make M101 ideal for 
studying the dependence of the aromatics on the local 
conditions, specifically metallicity and ionization. 

2. DATA 

The high sensitivity and spatial resolution provided by 
all three instruments on Spitzer allow individual HII re- 
gions in M101 to be studied in detail. Targeted spectro- 
scopic observations of a handful of HII regions and the 
nucleus were taken with the Infrared Spectrograph (IRS, 
iHouck et alj|2004bf ). Imaging of M101 was obtained with 
the Infrared Array Camera (IRAC, iFazio et aT1l2004h at 
3.6, 4.5, 5.8, 8.0 /jm and the Multiband Im aging Photom- 
etry for Spitzer (MIPS, iRieke et alj|200l at 24, 70, and 
160 /.tm. The IRAC and MIPS images allow the general 
behavior of the aromatics to be studied over the entire 
galaxy and the IRS spectra allow for the behavior of the 
aromatics to be probed in detail for selected HII regions. 
Past studies of M101 in the infrared have concentrated 
on radial gradients as the spatial resolution was not high 
eno ugh to study any but the few bright est HII regions 
(eg. Iffippelein etaTfll996t iPopescu et aT]|2005fh 

Preliminary resu lts from the IRAC and MIPS images 
were presented by IGordon et alj (l2006bf) and from the 
IRS spectroscopy in IGordon et al.l (|2006af) . This work 
focuses on the intense star forming HII regions in a sin- 
gle galaxy and is complementary to the Spitzer pro- 
gram on starburst galax i es led by C harles Engelbracht 
(jEngelbracht et alj 120051 120061 I2008D which focuses on 
the properties whole galaxies undergoing intensive star 
formation. 

2.1. IRAC and MIPS Images 

The IRAC images of M101 were obtained on 8 March 
2004. Two maps were taken separated by a few hours 
to allow for asteroid rejection. Each map consisted of 
13 x 13 tiles with 150"offsets (1/2 the array width) be- 
tween tile positions in both array dimensions . Each 
exposure was taken in high dynamic range (HDR) mode 
resulting in two images, one with a 10.4 sec exposure 
time and one with a 0.4 second exposure time. The im- 
ages were reduced with the Spitzer Science Center (SSC) 
pipeline S13.2.0 and combined using MOPEX program 
(v030106). As part of the mosaicking, overlapping re- 
gions of adjacent images were used to correct for bias 
drifts in each image. The calibration of the images was 
corrected to an infinite aperture (true surface brightness 
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Fig. 1. — IRAC, MIPS, and Ho images of M101 at full resolution. The IRAC 4.5 and 5.8 fim images have not been show n as they look 
very similar to the IRAC 3.6 and 8.0 fira images, respectively. The field-of-view of the images is 30' X 30'. The Ha image I) van Zee et al.l 
1998) is included to show the locations of HII regions. The streaking seen SE of M101 at 8 /im and roughly East- West at 70 and 160 (im 
are residual instrumental artifacts. 
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units) using the values given bv lReach et all (|2005t ) . The 
final mosaics have exposure times of ~85 seconds/pixel 
and are shown in Figure [T] 

The MIPS images of M101 were obtained on 10-11 May 
2004. Like the IRAC images, two maps were taken sep- 
arated by a few hours for asteroid rejection. Each map 
consisted of 14 medium rate 0.75 degree scan legs with 
a cross scan offset of 148". We use d version 3.06 of the 
MIPS Data Analysis Tool (DAT) ([Gordon et alj l2005h 
to do the basic processing and final mosaicking of the 
individual images. Extra processing steps before mo- 
saicking used programs written specifically to improve 
the reduction of large, well-resolved galaxies. At 24 /im 
the extra steps included readout offset correction, scan 
mirror dependent flat fields, array averaged background 
subtraction (using a low order polynomial fit to the data 
in each scan leg but excluding M101), and rejection of 
the 1st 5 images in each scan leg due to bias boost tran- 
sients. At 70 and 160 /im the extra steps were a pixel 
dependent background subtraction for each map (using 
a low order polynomial fit to the data in each scan leg, 
again excluding M101) and a spatial cosmic ray cleaning 
(160 /im only). The calibration of th e MIPS bands was 
done using the latest, official factors (|Engelbracht et alj 
[20071 iGordon et all [20071 IStansberrv et all 120071) . The 
5cr depth achieved in a single beam by these observa- 
tions is 0.050, 0.53, and 0.74 MJy sr" 1 for beam FWHM 
of 6", 18", and 40" and exposure times of -200, ~80, and 
~18 seconds/pixel for 24, 70, and 160 /im, respectively. 
The final mosaics are shown in Figure [1] The streaking 
seen along the scan direction at 70 & 160 /im is caused 
by residual instrumental artifacts and only produces a 
significant worsening of the noise at 70 /im. The on- 
line MIPS SENS-PET program predicts approximately 
the observed noise at 24 & 160 /im but a factor of 3 
lower noise than observed for the 70 /im observations. 
This difference is likely due to the residual streaking seen 
along the scan mirror direction which can be suppressed 
by combining scan map s taken at different scan angles 
([Tabatabaei et alJ l2007). Additional MIPS scan maps of 
M101 are planned and will be the subject of a future 
paper. 

2.2. IRS Spectra 

Low resolution IRS spectra from 5-38 /im of 7 HII re- 
gions and the nucleus were taken in the spectral map- 
ping mode with full slit width steps perpendicular to the 
slit. The locations of these 7 regions as well as other 
HII regions with measured electron tem peratures and 
metallicities from iKennicutt et al.l (|2003af ) are shown in 
Fig. [2] and their basic details are given in Table [T] The 
ShortLow spectra (5-15 /im) were taken with 6 steps 
and the LongLow (15-38 /im) were taken with 2 steps so 
that both spectral apertures covered the same regions. 
T he metallicities [12+lo g(Q/H)1 are taken from table 5 
of IKennicutt et all (|2003al ) except for the values for the 
nucleus and Hodge 602, which are the constant term in 
their equation 5. The ShortLow and LongLow obser- 
vations were taken approximately 3 months apart so the 
orientations of the two slits were approximately the same 
on the sky. 

These spectral mapping observations were reduced 
with the SSC pipeline (ve rsion S13/S14) and combined 
using CUBISM program ([Smith et all l2007al ). Spectra 



were extracted from the cubes for each HII region and 
the nucleus using an aperture with a radius of 10" and 
sky annuli from 12 to 18" referenced at 24 /im. The size 
of the aperture and sky annuli were varied linearly with 
wavelength to account roughly for the changing diffrac- 
tion limited point-spread-function (PSF), with a mini- 
mum object aperture radius of 5" to account of the finite 
size of the HII regions. The outlines of the individual IRS 
slits and extraction apertures are shown in Fig. 2] super- 
imposed on cutouts of the IRAC 8 /im and MIPS 24 /im 
images centered on each HII region. There are issues of 
crowding with nearby HII regions for some of the targets, 
but for the most part the targets are fairly well separated. 
The measured noise in the sky annuli region was used to 
compute the uncertainty on each point in the spectrum. 

The six spectral segments for each target were merged 
into a single spectrum using the overlap regions to adjust 
the level of each segment to be consistent. The overall 
level of the spectra was set by requiring each spectrum to 
be consistent with the MIPS 24 /im measurement given 
in Table [31 The main adjustment is at the transition 
between the ShortLow and LongLow spectra. This is 
not surprising given the observing strategy where the 
region is sampled with 6 ShortLow positions, but only 
2 LongLow positions. The resulting merged spectra are 
shown in Fig. [3J 

2.3. Resolution Matching 

One of the challenges of working with Spitzer data is 
the large change in PSF with wavelength as Spitzer is 
diffraction limited at wavelengths longer than ~6 /im. 
The FWHM (full width half maximum) of the Spitzer 
PSF varies from around 2" for the IRAC bands to 5", 
18", and 40" for the MIPS 24, 70, and 160 /im bands. To 
allow for accurate comparisons between wavelengths, it 
is important to convolve all images of interest to the low- 
est common resolution. Without this step, it is not easy 
to compare extended sources at different Spitzer bands 
given that the physical region probed changes with wave- 
length. The Spitzer PSFs are much more complex than 
simple Guassians with pronounced Airy rings and diffrac- 
tion spikes. This makes it definitely non-optimal to use 
simple Gaussian kernels to achieve the desired common 
resolution. It is also possible to put a set of images on a 
common resolution by convolving each image with all the 
other image PSFs, but this clearly creates a lower reso- 
lution common PSF than the method detailed below. 

The optimal method in the case of complex PSFs is to 
construct convolution kernels that transform the input 
PSF to the desired common resolution PSF. Theoreti- 
cally, the desired convolution kernel can be constructed 
using 2D Fast Fourier Transforms (FFTs) and is just 

K(x ,)=FFT- 1 ( FFT[PSF2(a: ^ )] l (1) 
\ FFT [PSFi (x, y)] J lij 

where PSFi is the input image PSF, PSF2 is the desired 
common resolution PSF, and K is the convolution ker- 
nel that transforms PSFi to PSF2. Unfortunately, the 
high frequency numeric noise in PSFi is greatly amplified 
and the resulting kernel is dominated by noise. The so- 
lution is to add a clamping function that attenuates the 
high frequency noise FFT [PSFi (a;, y)]. As the dominant 
signal we are interested in is radial, we adopt a radial 
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Hanning truncation function, Wh{ui), with the form 



18 fe) 



to < UJ Dl 

UJ > UJo 



(2) 



Fig. 2. — A 3 color image of M101 is shown annotated with the locations of HII regions with measured metallicities. The red, green, and 
blue channels are filled with the MIPS 24 fim, IRAC 8 fivci, and IRAC 3.6 /an images, respectively. The regions that were observed with 
IRS are identified as having a box around their names. 

This method of creating convolution kernels creates a 
set of images with a common PSF of the lowest resolu- 
tion in the set. For this method to work, the PSFs of 
each image have to be known to a reasonable degree and 
convolution kernels created for each unique PSF in the 
image set. The value of uj is set as large as possible 
for each kernel by visual inspection of the convolution 
of PSFi(a;, y) with K(x,y). The goal is for this con- 
volution to reproduce PSF2(a;,?/) with the highest pos- 
sible accuracy while minimizing the noise injected into 
the operation. An example showing the performance of 
the 24-to-70 /urn convolution kernel is given in Fig. O 
The close (but not perfect) correspondence between the 
70 /itm PSF and the convolution of the 24 jum PSF with 



where to is radial spatial frequency and to is the cutoff 
radial spatial frequency (|Brighaml [T988). The revised 
equation for the desired convolution kernel is then 



K(x,y) = FFT- 



FFT [PSF 2 (a;,y)] 
W A (w)FFT pSFifoy)] 



(3) 



The value of w Q in eq. [2] is chosen to be as large as pos- 
sible while still suppressing the high frequency noise in 
FFT [PSFi(x, y)]. 
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Fig. 3. — The IRS spectrum for all 8 regions with targeted IRS spectroscopy are shown. The left plot show the full spectra normalized 
to their average flux in the 20—30 (im region and on a log scale. The locations of aromatic and atomic emission lines are given. The right 
plot shows the 5-15 fim portion of the spectra normalized to their average flux at 10 /an and on a linear scale. The legend identifying the 
regions is given in this plot. 



TABLE 1 

M101 Regions with measured metallicities 



Name 


Other Names 




Coordinates 




R/R a 


12+log(0/H) 


IRS Spectra 


Nucleus 




14 


03 12.48 54 20 


55.4 


0.00 


8.76 ± 0.06 


X 


Hodge 602 




14 


03 10.22 54 20 


57.8 


0.02 


8.76 ± 0.06 


X 


Hodge 1013 




14 


03 31.39 54 21 


14.5 


0.19 


8.71 ± 0.05 




Searle 5 


Hodge 336 


14 


02 55.05 54 22 


26.6 


0.21 


8.55 ±0.16 


X 


NGC 5461 


Hodge 1105 


14 


03 41.36 54 19 


04.9 


0.32 


8.50 ± 0.03 


X 


NGC 5447 


Hodge 128 


14 


02 28.18 54 16 


26.3 


0.55 


8.45 ± 0.04 


X 


NGC 5462 


Hodge 1170 


14 


03 53.19 54 22 


06.3 


0.44 


8.37 ±0.07 


X 


NGC 5455 


Hodge 409 


14 


03 01.13 54 14 


28.7 


0.46 


8.25 ± 0.03 


X 


Hodge 67 




14 


02 19.92 54 19 


56.4 


0.54 


8.14 ±0.08 




Hodge 70/71 




14 


02 20.50 54 17 


46.0 


0.57 


8.10 ±0.12 




Searle 12 


Hodge 1216 


14 


04 11.11 54 25 


17.8 


0.67 


8.16 ±0.03 




NGC 5471 


NGC 5471-D 


14 


04 29.35 54 23 


46.4 


0.80 


8.09 ±0.03 


X 


Hodge 681 




14 


03 13.64 54 35 


43.0 


1.03 


7.92 ± 0.09 




a Ro = 14.42' 



the 24-to-70 /im kernel is a validation of the accuracy of 
this method. 

For the work on M101, we created convolution kernels 
that transform the IRAC PSFs to the MIPS 24 fmx PSF 
for use in §3.31 and kernels that transform the IRAC, 
MIPS 24 /im and MIPS 70 /im PSFs to the MIPS 160 fan 
PSF for use in H3A{ Ideally, we would have done the same 
to match the resolutions of each IRS wavelength to the 
resolution of the longest IRS wavelength. Unfortunately, 
the IRS PSFs are not nearly as well characterized as the 
IRAC and MIPS PSFs and our small regions mapped 
with IRS would introduce significant edge effects in any 
convolutions. Instead, we varied the extraction apertures 
to account for the changing IRS PSF with wavelength. 

3. RESULTS 

With the data presented in this paper, the behavior of 
the aromatic features in M101 can be probed both pho- 
tometrically and spectroscopically. The IRAC and MIPS 
total M101 fluxes are combined with previous measure- 
ments to probe the global SED of M101 f fl3~Tj) . Using 
the IRS spectra of the targeted HII regions, we probe the 
aromatic feature behavior in detail in £ 13.21 Returning to 
the imaging data, the behavior of the aromatics for an 
expanded sample of regions is studied ( £|3.3[) and, finally, 



the overall infrared morphology is used to gain further 
insight into the aromatics in M101 ( §3.4|) . Combining 
both the spectroscopic and photometric approaches pro- 
duces a richer picture of the aromatic feature behavior 
than is possible with either approach separately. 

3.1. Global SED 

The global IR SED of M101 is useful to put this galaxy 
in context with other galaxies and to probe the overall 
aromatic and dust content. The IRAC and MIPS global 
fluxes were measured using a circular aperture with a ra- 
dius of 13f3 and a sky annulus from 14.'2 to 16.'7. The 
uncertainties on these fluxes are dominated by the uncer- 
tainties in the absolute calibrations. The IRAC absolute 
calibration uncertainties are taken as 5% w hich is larger 
than the point source uncertainty of ~2% (|Reach et al.1 
l2005f) as the IRAC scattered light behavior makes the 
extended source calibration less accurate. The MIP S ab- 
solute uncertainties are 4% (|Engelbra£tolTaT][2llO^. 5% 
(jGordon et al.l[2007h . and 12% (|Stansberrv et alll200l . 
We have tabulated previous IR measurements of the to- 
tal M101 e mission in Table [21 We have included the 
ISOPH OT (jLemke et all [l99l) flux from iStickel et al l 
(|2004l) but not the fluxes given bv lTuffs k Gabriell (j2003f ) 
as these fluxes are known to be lower limits. The Mid- 
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NGC 5462 
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NGC 5455 





MIPS 34 M™ 



IRAC 8 Mm 



MIPS 24 M™ 

NGC 5447 
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IRAC 8 .iiiii 



MIPS 24 Mm 
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Fig. 4. — Cutouts centered on each of the 8 regions with targeted IRS spectroscopy are shown. Each cutout is 2' X 2' and shows the 
circular extraction aperture (and sky annulus) at that wavelength and the outline of the ShortLow and LongLow slits on the IRAC 8 /im 
and MIPS 24 fim images, respectively. For Hodge 602, the nearby nuclear region was masked (black square shown) to allow for an accurate 
sky measurement. 

simple stellar population and a three component dust 
grain emissi on model. The detail s of this modeling can 
be found in iMarleau et al.1 (|2006T ). While this model is 
not unique, it is illustrative of the dust properties im- 
plied by the SED. The global IR SED of M101 includes 
strong aromatic emission features as well as emission 
from warm (T - 56 K) and cold (T - 18 K) dust. 
The total dust mass of 1 x 10 s M Q is dominated by the 
coldest component and is much larger than the previ- 
ous estimate of 4.8 x 10 6 M based on IRAS data alone 
(jYoung et al.l fl989) . Using our new measurement of the 
dust mass, the M101 gas-to- dust ratio is ^300 using a 
HI mass of 2.5 x 10 10 M Q (|Bosma et al.lll981l) and a 



course Space Experiment (MSX, ICohen e t al. 200 l|) ob- 
serve d M101 in four mid-infrared bands (Kracmcr et al. 
2002), but only had a high signal-to- noise detection in 
MSX Band A. The MSX Band A flux was measured us- 
ing the same apertures as the IRAC and MIPS fluxes 
from the image available from NED. The uncertainty on 
this measurement was set to 20% to account for the dif- 
ficulty in sky subtraction for these data. 

The global IR SED of M101 is shown in Fig. © and 
looks s imilar to many o f the S pitzcr Nearby Galaxies 
burvey iKennicutt eFall (I2003U SINGS,) galaxy global 
SEDs pale et al.l [2005. 20071: iDraine et al.l 12007? ) . We 
fit the global SED of M101 with the combination of a 
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TABLE 2 
M101 Global Fluxes 





Wavelength 


Bandwidth 


Flux 






Band 


[Mm] 


H 


[Jy] 


Origin 


Reference 


J 


1.235 


0.162 


3.94 ± 0.12 


2MASS 


1 


H 


1.662 


0.251 


4 88 + 19 


2MASS 


1 


K 


2.159 


0.262 


4 i ft 4- n on 


2M ASS 


I 


IRAC1 


3.550 


0.681 


2.84 ± 0.14 


Spitzer /IRAC 


2 


IRAC2 


4.493 


0.872 


i 7fi + n f)9 


Snitzer/TR AC 


2 


IRAC3 


5.731 


1.250 


3.69 ± 0.18 


Spitzer /IRAC 


2 


LW2 


6.75 


3.5 


6 03 + 442 


ISO /TSOCAM 


3 


IRAC4 


7.872 


2.526 


7 2fi + (1 Sfi 


Snitzpr/TR AC 


2 


MSXA 


8.276 


3.362 


8.86 ± 1.77 


MSX 


2 


IRAS12 


12 


7.0 


6.2 ± 0.93 


IRAS 


1 


COBE12 


12 


6.48 


< 16.9 


COBE/DIRBE 


5 


LW3 


15.0 


6.0 


5.42 ± 0.501 


ISO/ISOCAM 


3 


MIPS24 


23.7 


4.7 


10.5 ±0.4 


Spitzer/MIPS 


2 


IRAS25 


25 


11.2 


11.8 ± 1.77 


IRAS 


i 


COBE25 


25 


8.60 


12.2 ±2.44 


COBE/DIRBE 


5 


IRAS60 


60 


32.5 


88.0 ± 13.2 


IRAS 


■1 


COBE60 


60 


27.84 


82.8 ±16.6 


COBE/DIRBE 


5 


MIPS70 


71.0 


19.0 


116.3 ± 7.0 


Spitzer/MIPS 


2 


IRAS100 


100 


31.5 


252.8 ± 37.9 


IRAS 


4 


COBE100 


100 


32.47 


259.0 ±51.8 


COBE/DIRBE 


5 


COBE140 


140 


39.53 


< 530.0 


COBE/DIRBE 


5 


MIPS160 


156 


35.0 


405.3 ± 35.0 


Spitzer/MIPS 


2 


ISO170 


170 


90.2 


495.8 ± 74.4 


ISO/ISOPHOT 


6 


COBE240 


240 


95.04 


< 497.5 


COBE/DIRBE 


5 



Ref ere nces. — (1) Uarrett et al.1 (I2003D ; (2 ) this wor k; (3) iRoussel et all 
(f200H ^: (4l lRice et all (|1988f ) f5) IOdenwald etaL1(|T99g ) f6) lStickel et al.li|2004ft 




1 10 1 00 

radius [arcsec] 



Fig. 5.— The radial profiles of the 24 /im PSF, the 24-to-70 fim 
kernel, the 70 (im PSF, and the convolution of the 24 fim PSF and 
the 24-to-70 fim kernel are shown. 

H2 mass of 5.1 x 10 9 M©. This new gas-to-dust ra- 
tio is m uch smaller than th e prev ious measurement of 
2440 by iDevereux fe Yound (| 1990ft . This is in accor- 
dance with their suspicion that there was an order of 
magnitude more dust which was too cold for IRAS to 
measure. The dust-to-gas ratio we measure is consistent 
with th at expected for a galaxy with a global metallicity 
of 8. 4 (|Kennicutt et al.l 2003at iMoustakas fc KennicuttJ 
[2006h from the work of IPraine et alj (j2007ft on SINGS 
global SEDs. 

The global obscured SFR can be measured using 
the total IR luminosity and the calibration given by 
iKennicutt et al.l (|2003aft . The total IR luminosity of 
M101 determined us ing the TIR (total IR) formulation 
()Dale fc Heloul 12002ft is 2.5 x 10 43 ergs s" 1 which gives 
a SFR of 1.1 Mq year" 1 . The global obscured SFR 



can also be measured using the MIPS 24 /im luminosity 
(XF(24 /im)47r<i 2 , where d = distance) which is 7.05 x 
10 43 ergs s" 1 . This gives a SFR of 1.8 and 1.1 Mp-, yea r" 1 
using the calibration s of lAlonso-Herrero et al.l (|2006ft or 
ICalzetti et~all (|2007ft . respectively. 



3.2. Spectroscopic Measures of Aromatics 

Measuring the strengths of the aromatic features in 
the IRS spectra is not straightforward given that many 
of the features are blended with each other as well as with 
atomic emission lines. Various methods have been used 
in the past that attempt to define continuum points be- 
tween features and measure consistently portions of the 
aromatic features, but these are prone to systematic er- 
rors as the level of blendin g changes wi t h the s trength of 
the features. Fortunately. ISmith et al.l (|2007bf ) provided 
a decomposition tool called PAHFIT to fit simultane- 
ously the aromatic features, atomic and molecular emis- 
sion lines, dust and stellar continuum, and the dust ex- 
tinction in IRS spectra. Some previous w ork in this area 
has u sed similar fitting techniques (e.g., iMadden et al.l 
12006ft . but PAHFIT is the first such tool to fit simul- 
taneously all the components of the mid-IR spectra of 
star forming region. As PAHFIT takes into account the 
blending of the aromatic and atomic features, it produces 
fairly precise measurements of their individual strengths. 
We used the PAHFIT program to fit all of our M101 
spectra and the resulting fits for two of the regions are 
shown in Fig. [7J The atomic emission line strengths from 
the fits are given in Table [H The aromatic emission line 
strengths and equivalent widths are given in Tables[5]and 
El respectively. The line strengths are computed by in- 
tegrating the fit profile for the specific line and does not 
include any contribution to the continuum. The equiva- 
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Fig. 6. — The IR SED of M101 is shown from the new IRAC and MIPS observations presented here as well as measurements taken 
from the literature. The COBE/DIRBE upper limits are 3cr. The horizontal error bars give the widths of the bands. A simple stellar 
population plus 3 component dust emission model fit is shown. Each dust component is modeled with a single grain size and the temperature 
given is the equilibrium temperature. The PAH component is illustrative only as it is assumed to be in thermal equilibrium which is an 
approximation as this component is known to be heated stochastically. 



lent width is 

where I(v) = I{v) line + I(u) con t (both from the PAHFIT 
results) and is defined such that the equivalent width is 
positive for an emission line. The uncertainties quoted 
are returned from PAHFIT and are the formal uncertain- 
ties on the quantities. These PAHFIT returned uncer- 
tainties are based on the input spectrum uncertainties 
which in turn arc defined using the background noise. 

There are a number of diagnostics of the physical con- 
ditions in the HII regions contained in the IRS spectra 
themselves as well as literature optical spectral stud- 
ies. The IRS spectra give probes of the ionization 
([NeIII]/[NeII], [SIV]/[SIII], [ Arlll]/ [Aril]) and e lectron 
density ([SIII] 18/[SIII] 33) llGiveon et al.ll200l . The 
ionization probe is related to the radiation field hard- 
ness, but it is not easy to make a direct translation. 
The optical studies provide high quality metallicities 
[log(0/H)+12], electron tempe ratures, and ionization 
(e.g., [01111 5 007/[OII1 3727). (IKennicutt et all l2003al : 
iBresoliij 12006? ). Finally, the MIPS 24 /im flux provides 
a rough probe of the radiation field density conv olved 
with the amount of dust present. IWu et all (|2006f ) used 
a similar measurement, L22[im/V, where L22fj.n1 was de- 
termined from IRS Red Peakup images and V is the vol- 
ume of the region probed (in this case the whole galaxy) . 
This was described as a probe of the UV radiation field 
density assuming dust absorbs most of the energy and 
there is a fairly constant conversion from L22/im and to- 
tal IR flux. In our case, the V is similar for all our 
HII regions because we are probing regions in the same 
galaxy with photometry performed with the same sized 
apertures. 



The different diagnostics are not all independent and 
Fig. [5] illustrates this point. The correlations among the 
different ioniza tion probes were see n in ISO observations 
of HII regions dGiveon et alll2002[ ) and starburst galax- 
ies (|Verma et al.l l2003f ) . Not surprisingly, we find the 
same for M101 HII regions for the correlation between 
[NeIII]/[NeII] and [SIV]/[SIII] 18 (Fig. pi, [ArIII]/[ArII] 
is not plotted as both lines are only significantly detected 
in one of our HII regions). The infrared ([NeIII]/[NcII]) 
and optical (dereddened [OIII] 5007/ [Oil 3727]) probes 
of ionization are seen to be correlated (Fig. [SJd), but 
not nearly as well as the two infrared ionization probes 
(Fig. [5^). Metallicity and ionization are known to be 
roughly correlated. This rough correlation is caused by 
decreasing line blanketing and harder spectra in hot stars 
as the metallicity decreases. In the M101 HII regions, we 
also see these two quantities are correlated (Fig. [8k), but 
the correlation has a larger scatter than can be accounted 
for by the observational uncertainties. Fig. [Hp 1 illustrates 
that all the HII regions have similar electron densities, 
< 100 cm " 3 , as determined fr om the [SIII] 18/ [SIII] 33 
line ratio dGiveon et "ail |2002[) a s well as optical spec- 
troscopy (jKennicutt et al.ll2003al) . This is likely due to 
the different ionization energies between [Oil] (13.6 eV) 
and [Nell] (21.6 eV) and [SIII] (23.3 eV). The rough cor- 
relation between the optical and infrared ionization mea- 
sures indicates that there may well be significant varia- 
tion in the radiation field hardness among the different 
HII regions. The radiation field intensity as probed by 
the MIPS 24 flux is not correlated with any of the 
physical quantities (Fig. [8^) although this lack of cor- 
relation is mainly based on the large 24 fim flux seen 
for NGC 5461. Finally, the optically determined elec- 
tron temperatures are well correlated with metallicity 
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Fig. 7. — The IRS spectra for the Nucleus and NGC 5447 are shown along with the PAHFIT results. The total PAHFIT spectrum is 
shown as well as the decomposition of the total (red) into dust continuum (blue), atomic emission lines (purple), and aromatic features 
(green). The silicate absorption has been applied to the PAHFIT results. No significant starlight is found by the PAHFIT program. The 
dust continuum is fit with up to 8 com ponents with fixed temperatures and is meant to produce a continuum which is physical motivated 
(sec section 4.13 of Smith ct al. 2007b). This continuum fit is not unique and the strengths of the different components do not necessarily 
correspond to different real grain components. 



(Fig. [8f). Overall, there seem to be three quantities 
that probe different physical properties of the HII re- 
gions. They are log(0/H)+12 (metallicity and electron 
temperature), [NeIII]/[NeII] and [SIV] / [SIII] 18 (ioniza- 
tion), and MIPS 24 /im flux (radiation field intensity). 

Given that the [NeIII]/[NeII] and [SIV] /[SIII] 18 
ratios are well correlated, it is possible to combine 
the two ratios into a composite ionization measure- 
ment with smaller unc ertainties. In this paper and 
lEngelbracht et alJ (2008), we propose a composite mea- 
sure called the ionization index (or II for short). It is de- 
fined as the weighted average of log([SIV]/S[III] 18) con- 
verted to the [NeIII]/[NeII] scale and log([NeIII]/[NeII]). 
The conversion of [SIV]/S[III] 18 to the [NeIII]/[NeII] 
scale is done using the observed correlation for the M101 
HII regions and t he st arburst galaxies presented by 
lEngelbracht et all (|2008h . This correlation for just the 
M101 HII regions (Fig. [Sk.) is 

log (y) = (0.93 ± 0.11) + (1.85 ± 0.34) log (x) (5) 

where x = [SIV] /[SIII] 18 and y = [NeIII]/[NeII]. The 
combined correlation of the M101 HII regions and the 
starburst galaxies is 



log (y) = (0.71 ± 0.08) + (1.58 ± 0.13) log (x) 



(G) 



We adopt the combined correlation to determine the II 
values for the M101 HII regions. Thus, II can be calcu- 
lated from 



(if 



0.71 + 1.58 log 



[SIV] 
[SIII] 18 



(7) 



We use the II values for the rest of this paper instead of 
the [NeIII]/[NeII] or [SIV] /[SIII] 18 ratios as the II has 
less noise and is available for all the M101 HII regions. 

We have chosen to call the combined measurement of 
[NeIII]/[NeII] and [SIV] /[SIII] an ionization index as this 
is the physical measurement which is being made. Often 
these two ratios are used as an indicator of the radia- 
tion field hardness (or T e ff). These ratios are dependent 
on the metallicity, ionization parameter and morphol- 
ogy, and age of an HII region. The dependence is such 
that the ratios decrease as the metallicity increases, the 



ionization parameter decreases, and the age decreases 
(jRigby fc Riekd I2004T ). This explains why these ratios 
are only roughly correlated with metallicity. The scatter 
in the metallicity versus [Nelll] / [Nell] ratio (Fig. [St) is 
likely due to different HII region ages and/or ionization 
parameters. Disentangling the different effects is compli- 
cated and beyond the scope of this paper as we are only 
interested in these ratios as a measure of the amount 
of processing which might have taken place in the HII 
regions. 

The equivalent widths of the aromatic features are 
plotted versus metallicity and II in Fig. [9l The equiv- 
alent width measures the strength of the aromatic fea- 
tures versus the underlying mid-infrared dust continuum 
emission. As both the aromatic features and underly- 
ing continuum are measured at the same wavelengths, 
they are likely to have undergone similar excitation and, 
therefore, the equivalent width is a measurement of the 
abundance ratio of the aromatic carriers to small dust 
grains. Visually, it appears that the equivalent width of 
the aromatic features is better correlated with II than 
with metallicity. We have quantified this by fitting the 
measurements with two different functions. The first is 
a simple power law 

y = ax b (8) 

where x = 10 n and y is the normalized aromatic equiv- 
alent width. The second equation is a power law plus a 
constant with the form 

y=[ax b + c]- 1 . (9) 

The best fits using both of these equations are shown in 
Fig. M The reduced x 2 values for the data versus metal- 
licity are 356 and 339 for the simple power and power 
law plus a constant, respectively. The reduced x 2 val- 
ues for the data versus II are 240 and 195 for the simple 
power and power law plus a constant, respectively. This 
indicates that the aromatic feature equivalent widths are 
better correlated with ionization than metallicity. In ad- 
dition, the functional form favored is the power law plus 
constant (Eq. [9]). The behavior of the aromatic equiva- 
lent widths can be described as constant up to a thresh- 
old II (~ which is a [Nelll] / [Nell] value of 1) and then 



Aromatics in M101 



11 



10.0 



1.0 



: a) 








' 7^ 







0.1 



1.0 



[SIV] / [SHI] 18 



8.8 



8.6 



B.'l 



8.2 

















- c) 





0.1 



1.0 

[Nell I] / [Nell 



10.0 



10.0 



1.0 



e) 



100 

F„(24 /im) [mjy] 



1000 



0.35 [ 



0.30 ■ 



b) 



rn 0.25 - 



0.20 • 



0.15 



0.10 [ 



0. « 



1.0 10.0 
[Nelll] / [Nell] 



d) 



10.0 r 



0.1 r 



0.1 



' .0 



[Olll] 5007 / [Oil] 3727 



8.4 - 



8.2 - 



f) 



7.0x10 3 8.0x10 3 9.0x10 3 1.0x10 4 1.1 x10* 1.2x10 4 
T e (0 + ,N + ,S + ) 



Fig. 8. — Diagnostics of the physical conditions in HII regions are plotted versus each other to illustrate their interdependencies. These 
probes are the ionization ( [Nelll] / [Nell] and [SIV] / [SIII] 18, dereddened [OIII] 5007/ [Oil] 3727), electron density ([SIII] 18/[SIII] 33), 
metallicity [log(0/H) + 12], the optical emission line electron temperature [T e (0 + ,N + ,S + )], and radiation field intensity (MIPS 24 fim 
flux). 



a reduction in equivalent width with a power law form 
at higher II values. The best fit between the II and nor- 
malized equivalent width measurements is 

y = [0.55a; - 98 + 0.73] _1 . (10) 
The expected values for the equivalent widths for a par- 
ticular aromatic feature can be determined from this 
equation multiplied by the average equivalent widths 
that are given in the caption of Fig. [5] 

The ratio of the integrated strengths of different aro- 
matic features is a diagnostic of the physical mechanism 
of processing. For example, models using PAHs as the 
carriers of the aromatic features predict the strength of 



the 6.22, 7.70, and 8.61 /jm features to strengthen rel- 
ative to the 11.3 and 12.7 /im features as the ioniza- 
tion level inc reases or to weaken as the smaller P AHs 
are removed (jBakes et alj|2001t [Peeters et all 120021) . In 
Fig. [TO] we show the ratio of the 4 strongest aromatic 
features to the 11.3 ^m feature versus metallicity and II. 
The ratios do not vary significantly from object to ob- 
ject nor does there seem to be a significant trend with 
either of the two diagnostics. While there is a large scat- 
ter in the ratios of around a factor of two, the large 
measurement uncertainties mean that the ratios could 
be consistent with no change with either of the two di- 
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Fig. 9. — The normalized equivalent widths of the 5 highest equivalent width aromatic features (AFs) are plotted versus metallicity (left) 
and ionization index (II [see £|3. 2 1 11, right). The normalization was done to the average equivalent width of each AF. The normalization 
values are 0.87, 3.01, 1.04, 1.32, and 0.56 fim for the 6.22, 7.70, 8.61, 11.30, and 12.70 f im a romatic features, respectively. The dotted and 
dashed lines give the best fit power law and power law plus constant, respectively (see £|3.2| I, 
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Fig. 10. — The ratio of the strength of the 4 strongest aromatic features (AFs) to the 11.3 AF are plotted versus metallicity (left) and 
ionization index (II, right). The normalization was done to the average ratio of each AF. The normalization values are 1.01, 3.41, 0.80, and 
0.57 for the 6.22, 7.70, 8.61, and 12.70 /an aromatic features, respectively. 



agnostics. This is different than the preliminary result 
reported bv lGordon et alj (|2006af ) where significant vari- 
ation in the ratios was found. We have improved both 
the IRS data reductions and aromatic feature equivalent 
width measurements, which accounts for the differences. 
The lack of ratio variations is interesting in light of the 
strong weakening of the aromatics above a threshold in 
II (Fig. [5]). The aromatic features weaken by over a fac- 
tor of 10, but the ratio of different aromatics does not 
change by more than a factor of 2. 

3.3. Photometric Measures of Aromatics 

While IRS spectroscopy is clearly the most accurate 
way to study the aromatic features, it would be useful 
to be able to probe the aromatic feature strengths using 
some combination of IRAC and MIPS imaging. Having 
a photometric probe to the aromatic features would al- 
low for a larger range of objects and environments to 
be studied given that imaging is faster to acquire and 
can probe regions that are too faint for spectroscopy. 
lEngelbracht et all (|2005f ) proposed the combination of 
8-to-4.5 and 8-to-24 /xm flux density ratios as just such 
a measure of the aromatic feature strength given that 
the IRAC 8 yum band probes the strong 8 /jm aromatic 



complex (composed of the 7.7, 8.3, and 8.6 /xm aromatic 
features) and the IRAC 4.5 /im and MIPS 24 /im bands 
probe different parts of the small grain dust emission. 

The 8-to-4.5 and 8-to-24 flux density ratios are 
measuring something very similar to an aromatic fea- 
ture equivalent width as it can be thought of as the ra- 
tio of the aromatic emission divided by the underlying 
small grain emission. It is not the same as an equivalent 
width as the small grain emission measurement is done 
at shorter and longer wavelengths (4.5 and 24 /xm, re- 
spectively, instead of 8 /im) than is done when using the 
IRS spectroscopy. In other words, the 4.5 /im flux probes 
hotter and 24 /xm cooler dust grains than the 8 /im spec- 
troscopic continuum. Thus, variations in 8-to-4.5 and 
8-to-24 /im can be due to changes in the dust emission 
continuum instead of the 8 /im aromatic complex. Yet 
if both ratios change in a consistent manner, it is highly 
likely it is the 8 /xm aromatic complex which is varying. 
In Figure [TT] we plot the 8-to-24 /im versus 8-to-4.5 /im 
flux density ratios for a sampl e of M101 HII regi o ns an d 
a sample of starburst galaxies lEngelbrac ht et all ([2008) . 
The M101 HII region photometry was performed on the 
MIPS 24 /xm and IRAC 3.6, 4.5, & 8.0 images convolved 
to the 24 /xm PSF ( q2.3[) . The photometry was mea- 
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Fig. 11. — The f„ s (8 /im)//^ s (4.5 /im) ratio is plotted versus 
the f^ s (8 /im) / J '~ :s (24 fim) ratio where the superscript "SS" des- 
ignates these fluxes as stellar subtracted. The contributions to the 
measured fluxes from stellar emission was subtracted by scaling 
the measured 3.6 (im flux by 0.57, 0.23, and 0.031 for the 4.5 fira, 
8.0 /im, and 24 (im bands, respectively. The HII regions are plotted 
as so lid, red circles and the starburst galaxies CEiiKclbrac bt et al.l 
2008) as open, green squares. The error bars are derived from the 
measurement uncertainties and include the contribution from stel- 
lar subtraction. The dotted line shows a spectroscopically-chosen 
color separation between galaxies with (o utside the box) and with - 
out (inside the box) aromatic emission (Engclbracht et al. 2008). 

sured using a 10" radius object aperture and a back- 
ground annulus with radii from 13 to 20"and is given 
in Table [3] The IRAC 3.6 /im band is dominated by 
stellar emission and was used to subtract the stellar 
contribution from the lo nger wavelength measurements 
(lEngelbracht et all 12005). Th e starburst photometry is 
from lEngelbracht et al.l ( 2008f ) and was also corrected for 
stellar emission. 

The plot shown in F ig. fTTI is similar to Fig. 1 of 
lEngelbracht et all (|2005h except that we have plotted 
the 8-to-4.5 /zm instead of 4.5-to-8 /im flux ratio so that 
both the ratios plotted behave in a similar sense where 
larger ratios imply stronger aromatic features. In ad- 
dition, we have subtracted the stellar contribution from 
all three bands involved using scaled 3.6 /im flux densi- 
ties. These two ratios probe the behavior of the 8 /tm 
aromatic complex plus underlying dust continuum emis- 
sion versus shorter wavelength dust continuum emission 
(8.0/4.5) and longer wavelength dust continuum emis- 
sion (8/24). Fig. mi shows that the M101 HII regions 
and starburst galaxies show the same behavior and cover 
similar ranges in both ratios. This is not surprising given 
that the simplest model of a starburst is a collection of 
HII regions. However, it is important to confirm that 
HII regions have mid-IR SEDs similar to those of star- 
burst galaxies and, as a result, are reasonable analogs for 
studying dust emission from starburst galaxies. 

Given that we have spectroscopic and photometric 
measurements of the 8 /jm aromatic feature complex, 
we can see how the photometric measure compares with 
the more accurate spectroscopic measurement. Figure li"2l 
gives this comparison where the measured spectroscopic 
equivalent width of the 8 /im aromatic complex is plot- 
ted versus the 8-to-4.5 /im and 8-to-24 /zm flux density 
ratios. It is clear that these two flux ratios are rough 
measures of the 8 /im complex. The 8-to-24 /tm ratio 



is probably a better measure given that the relationship 
is more linear than that using the 8-to-4.5 /im ratio and 
the measurement uncertainties are smaller. 

It is possible to combine both ratios in a way to create 
a better measure of the 8 /im complex. The basic idea 
is to use the 4.5 and 24 /zm measurements to predict 
the continuum at 8 /im due to small grain emission and 
combine this continuum prediction with the measured 
8 /tm flux to create a photo metric measure of th e 8 /im 
equivalent width. Following lEngelbrac ht et ail (|2008h . 
wc find that it is best predicted using 



F !/ (8/xm) cont - F„(4.5/zm) - 66 x F„(24/mi) 



0.34 



(11) 



which results from assuming the small grain emissio n 
follows a power law. Unlike Engclbrac ht et ail (2008), 
we have not used the stellar subtracted fluxes to deter- 
mine ^(8 /im) cont as we do not have K band fluxes to 
help measure the stellar flux and the stellar contribution 
is expected to be small for these resolved HII regions. 
Using the stellar subtracted fluxes produces systemati- 
cally larger photometric than the measured spectroscopic 
equivalent widths. The photometric equivalent width is 
then 



EQW(8 /im) p i 

1( 



F„(8/mi) 



1 I AA. (12) 



where 5X = 2.69 /im and is the width of the IRAC 8 /im 
band. This assumes that the 8 /tm complex fills the IRAC 
8 /tm band which is the case for many of the HII regions, 
but not all of them (see Fig. [3]). We have computed the 
photometry equivalent widths for the M101 regions and 
plotted them versus the spectroscopic equivalent widths 
in Fig. 1131 The photometric equivalent width measure- 
ment produces a good rough measurement of the 8 /tm 
aromatic co mplex for the M101 region s. A similar result 
is found by lEngelbrac ht et all (|2008h for the starburst 
galaxy sample. 

3.4. Morphology of Aromatic Emission 

The IRAC and MIPS images of M101 can be used to 
probe the morphology of M101 in the 8 /tm aromatic 
complex by applying the photometric measurement tech- 
niques introduced in H3.3I Fig. Q3] shows the simple 8-to- 
24 /im flux density ratio image and the EQW(8 /im)ph oto 
image. The IRAC images were convolved to the 24 /im 
resolution (i j2.3p of 6" (0.2 kpc at 6.7 Mpc) before com- 
puting either composite image. Caution should be exer- 
cised in interpreting these images as both involve stellar 
subtracted images. The stellar subtraction is done using 
a multiplicative scaling of the IRAC 3.6 /im image and 
the accuracy of this method has only been shown for in- 
tegrated fluxes. This method will oversubtract regions 
where the emission in the 3.6 /im image is dominated 
by dust emission instead of stellar emission. While the 
EQW(8 /im) p hoto image is expected to be a more accu- 
rate estimate of the true 8 /tm aromatic complex equiv- 
alent width than the 8-to-24 /im image, it suffers from 
higher noise given the sensitivity of this measurement to 
the more uncertain stellar subtracted 4.5 /im flux and 
8 /tm instrumental residuals. Both images show that the 
aromatic features are depressed in HII regions (white re- 
gions in Fig. [I}. This is illustrated in Fig. [T5] where 
the Ha image is shown in blue and the 8-to-24 /im and 
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TABLE 3 
Photometry for M101 regions 



name 


IRAC1 


IRAC2 


IRAC3 


IRAC4 


MIPS24 




3.6 fira 


4.5 fim 


5.8 lira 


8.0 fiar 


24 /im 




[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 



Nucleus 


15.00 


± 


0.119 


9.70 


± 


0.083 


19.20 


± 


0.222 


48.90 


± 


0.675 


96.20 


± 


1.221 


Hodge 602 


3.09 


± 


0.595 


1.94 


± 


0.364 


3.78 


± 


0.662 


9.35 


± 


1.545 


23.20 


± 


3.654 


Hodge 1013 


1.59 


± 


0.051 


1.47 


± 


0.035 


7.33 


± 


0.166 


20.00 


± 


0.447 


112.00 


± 


0.531 


Searle 5 


1.17 


± 


0.017 


1.00 


± 


0.012 


4.82 


± 


0.039 


13.20 


± 


0.126 


36.20 


± 


0.128 


NGC 5461 


7.01 


± 


0.143 


6.96 


± 


0.121 


28.40 


± 


0.534 


82.00 


± 


1.392 


923.00 


± 


5.994 


NGC 5462 


1.37 


± 


0.047 


1.16 


± 


0.040 


4.08 


± 


0.086 


10.40 


± 


0.213 


74.50 


± 


0.884 


NGC 5455 


1.92 


± 


0.009 


1.75 


± 


0.009 


5.36 


± 


0.031 


13.10 


± 


0.082 


134.00 


± 


0.388 


NGC 5447 


1.75 


± 


0.104 


1.47 


± 


0.091 


4.62 


± 


0.342 


12.30 


± 


0.898 


106.00 


± 


6.463 


Hodge 67 


0.19 


± 


0.008 


0.17 


± 


0.006 


0.43 


± 


0.012 


1.02 


± 


0.026 


6.94 


± 


0.071 


Hodge 70/71 


0.26 


± 


0.008 


0.26 


± 


0.006 


0.52 


± 


0.017 


1.23 


± 


0.043 


10.20 


± 


0.120 


Searle 12 


0.45 


± 


0.007 


0.44 


± 


0.005 


0.79 


± 


0.009 


1.76 


± 


0.019 


16.30 


± 


0.079 


NGC 5471 


1.56 


± 


0.012 


1.80 


± 


0.011 


2.10 


± 


0.019 


4.66 


± 


0.035 


122.00 


± 


0.777 


Hodge 681 


0.12 


± 


0.008 


0.10 


± 


0.005 


0.11 


± 


0.005 


0.14 


± 


0.004 


1.53 


± 


0.042 
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Fig. 12.— The F^ s (8 lira) / F^ s (4.5 fim) and F^ s (8 iiar)/F^ s (24 lira) ratios are shown plotted versus measured 8 lira aromatic feature 
complex equivalent width. 
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Fig. 13. — The photometric (sec i|3,3H versus spectroscopic mea- 
surements of the 8 fira aromatic complex equivalent width are plot- 
ted. The dashed line shows the one-to-one relationship. 

EQW(8 /im)photo images are shown in yellow. The pre- 
ponderance of yellow and blue in the first image shows 
that the 8-to-24 /im is anti-correlated with Ha emission. 
This is an indication that the 8 /im flux is depressed 
and/or the 24 /mi flux is enhanced in HII regions. The 
EQW(8 /im) p hoto+Ha image directly probes the strength 



of the 8 /xm aromatic complex versus Ha and shows that 
the separation between the yellow and blue is not as 
clean as in the 8-to-24 /im+Ha image. Thus, the en- 
hancement of 24 /im emission in HII regions is a signif- 
icant contributer to the morphology of the 8-to-24 /im 
image. In other words, the dust in HII regions is hot- 
ter (emitting more at 24 /im) than the surrounding dust, 
but this hotter dust does not have a corresponding in- 
crease in the 8 /im continuum or aromatic emission to 
keep the 8-to-24 /im ratio constant. But there is clearly 
also a change in the 8 /im emission, given that there are 
large HII complexes (e.g., NGC 5462, NGC 5447) in the 
second image which show the expected morphology with 
aromatic emission (yellow) surrounding the Ha emission. 

Finally, the aromatic feature and underlying contin- 
uum emission can be compared to that of the large grains 
by including the longer wavelength MIPS images in the 
analysis. Fig. [TBI shows two images: the first is the ratio 
of the 8 /im flux to the total infrared (XF(8 /im) /F{IR)) 
and the second is the ratio of the 24 /im flux to the total 
infrared (AF(24 /im) / F '(I R)) . Generally, the 8 /im emis- 
sion traces the aromatic features (with a minor contri- 
bution from the continuum), the 24 /im emission traces 
the hot dust emission from small grains, and the total 
infrared traces the bulk of the dust grains which dom- 
inate the dust mass and have a peak emission around 
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Fig. 14. — The 8-to-24 and EQW(8 ^tm) p hoto images of M101 are shown. The images are displayed in reverse B/W using a asinh stretch 
with ranges from 0.1-2.0 and 3-30 and pivots of 0.5 and 5.0, respectively. The inner 20' of M101 is shown. Only pixels with measurements 
above the sky noise are displayed. The EQW(8 /im) plloto image shows significant noise due to instrumental signatures present in the IRAC 
8 /im mosaic (see Fig. [TJ. 




Fig. 15. — The 8-to-24 and EQW(8 fJ.m) p ^ oto images of M101 are shown in yellow along with the Ho image in blue. The inner 20' of 
M101 is shown. Only pixels with measurements above the sky noise are displayed. 



160 /jm (see Fig. [BJ. To first order, XF(8 fim)/F(IR) 
and XF(24 fim)/F(IR) trace the ratio of aromatic grain 
mass to dust mass and small grain mass to dust mass, 
respectively, modified by the illuminating radiation field. 
All the images used were first convolved to the 160 /im 
resolution of 40" (1.3 kpc at 6.7 Mpc) using the kernels 
described in §2.31 The total infrared flux was calculated 
by a simple integration of the IRAC 8 /im through MIPS 
160 /Ltm images. The AF(8 fxm)/F(IR) image is domi- 
nated by a high ratio in the inner region with a fairly 
abrupt drop to much lower values. The brightest HII re- 
gions are seen to have shallow depressions compared to 
other regions at the same radii. The XF(24 /im) /F(IR) 
image displays a different morphology with a relatively 
flat ratio overall with strong enhancements in the HII re- 



gions. The depression of the aromatic emission is not due 
to lack of exciting photons as these are clearly present 
given the enhancement of the 24 /jm emission. The de- 
pression of the aromatic emission in HII regions must 
then be due to processing of the carriers of these fea- 
tures. The behavior of these two ratios supports the 
results of the comparison to the Ha given in the previ- 
ous paragraph. The variations of the 8 /im emission are 
dominated by variations in ionization while the 24 /im 
emissions are dominated by the radiation field intensity 
(as traced by the Ha emission). 

4. DISCUSSION 

The behavior of the aromatic feature equivalent widths 
versus metallicity and II seen in this paper for M101 
HII regions is also seen for starburst galaxies. This is 
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Fig. 17. — The normalized equivalent widths of the 5 highest equivalent width aromatic features (AFs) are plotted versus metallicity 
(left) and ionization ind ex (II, right). These two plots are the same as shown in Fig. [9] but now include the starburst sample from 
lEngelbracht et al . (2008) i n grey. The loca tions of the two most metal poor galaxies (SBS0335-052E and IZwl8) are labeled. In addition, 
the II ranges probed in lBrandl et aLl <f2006h and lWu et all ([2006 ) are indicated. The normalization was done the same as in Fig. [9] 



shown in Fig. \T7\ which combines the measurements for 
the M101 HII regions presented in t his paper with the 
starbu rst galaxy measurements from lEngelbracht et all 
(2008). Both samples show a better correlation of the 
aromatic equivalent widths versus II than metallicity. 
The combined M101 HII region and starburst galaxy 
samples probe ionization indexes between -1.5 and 1.4 
(equivalent to [NeIII]/[NeII] ratios of 0.03-25) and metal- 
licities (log(0/H)+12) from 7.9-8.8. There are two 
galaxies (SBS0335-052E and IZwl8) which have metal- 
licities below 7.9, but this is not enough points to statis- 
tically probe the behavior in this range. It is interesting 
to note that these two galaxies do not have the most ex- 
treme II values and fit fairly well in the trend of II versus 
aromatic feature equivalent widths. The residuals of fit 
of II versus aromatic feature equivalent width were plot- 
ted versus metallicity to see if the scatter around the fit 
was due to metallicity, but no trend was present. 

The correlation of the aromatic equivalent widths with 
II is not a simple power law, but shows a change at 



around a II value of zero. Below this value, the equiva- 
lent widths are fairly constant indicating that the mate- 
rial responsible for the aromatic features is as robust as 
the material responsible for the underlying dust contin- 
uum emission. Above this value, the equivalent widths 
drop significantly with increasing II. This implies that 
the material responsible for the aromatic features is be- 
ing modified/destroyed at a faster rate than the material 
responsible for the underlying dust continuum emission. 

The complex behavior of the aromatic equivalent 
widths with II explains the seemin gly contradictory re - 
sults found by previous authors. iBrandl et al.l (|2006l ) 
studied low i onizations ( —1.25 < II < 0.05) and found 
no variation. IWu et all (|2006f) studied high ionizations 
(—0.3 < II < 1.3) and found large variations. Only by 
studying the full range of II values does the complex 
nature of the correlation between II and aromatic equiv- 
alent width emerge. 

The complex behavior of the aromatic equivalent 
widths also explains why the bright point sources (e.g., 
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Fig. 18. — The ratio of the strength of the 4 strongest aromatic 
features (AFs) to the 11.3 AF are plotted versus the 6.2/11.3 AFs. 
The normalization was done to the average ratio of each AF. The 
normalization values are 1.01, 3.41, 0.80, and 0.57 for the 6.22, 7.70, 
8.61, and 12.70 £tm aromatic features, respectively. The points for 
the 6.22 AF defined the one-to-one correspondence. 



HII regions) in M101 change fairly abruptly from yellow 
to red at approximately the same radii in Fig. O The col- 
ors of the HII regions are primarily composed of green 
(IRAC 8 /mi) and red (MIPS 24 /mi). In the inner re- 
gions of M101, the HII regions have high metallicity and 
low ionization (e.g., Fig. [8k). This means that they have 
II values below and, thus, normal aromatic equivalent 
widths (e.g., IRAC 8 /im emission). In the outer regions, 
the HII regions have low metallicity and high ionization. 
Thus, they have much weaker aromatic emission which 
strongly reduces the green color. Red color in this partic- 
ular type of 3 color image is a strong indication of weak 
aromatic emission. 

Our finding that the aromatic equivalent widths in 
the M101 HII regions and starburst galaxies correlates 
better with a measure of processing (ionization) than 
with a measure of formation (metallicity) is not too sur- 
prising. There has been ample evidence for such pro- 
cessing in spatially resolved studie s in reflection nebulae 
and HII regions in the Milky Way (|Cesarskv et alj|200d: 
iBerne et al.l l2007h . But until this work, it is not been 
clear if the spatially resolved processing of the aromatic 
carriers also applies to the global spectra of massive star 
forming regions (i.e., HII regions and starburst galaxies). 
This work indicates that the processing of the aromatic 
carriers is happening on a global scale in massive star 
forming regions. 

Models where PAH molecules are the carriers of the 
aromatic features would predict a correlation of the ratios 
of different aromatic features with ionization index if the 
processing mechanisms were ionization or size selective 
destruction. The lack of such a correlation (see Fig. IT0|) 
seems to imply that these two processing mechanisms 
may not be the dominant cause of the strong correlation 
seen with ionization index. Fig. [18] plots the ratio of the 
strongest 4 aromatic features and the 11.3 /im aromatic 
versus the 6.2/11.3 aromatic ratio . In PAH models of the 
aromatic features ([Galliano! l2006h . the 6.2 /im aromatic 
traces small, ionized PAHs and the 11.3 /im aromatic 
traces larger, neutral PAHs. Thus, the 6.2/11.3 ratio 
traces ionized/neutral and/or small/large PAHs. The 



7.70 /im and 8.61 /tm aromatics also trace small, ion- 
ized PAHs and the 12.7 /im larger, neutral PAHs. While 
there is scatter in Fig. QH the 7.70/11.3 and 8.61/11.3 
ratios correlate fairly well with the 6.2/11.3 ratio. On 
the other hand, the 12.7/11.3 does not correlate at all 
with the 6.2/11.2 ratio. Thus, the behavior of the rel- 
ative strengths of the 5 strongest aromatic features is 
consistent with the PAH molecule model. The lack of 
correlation seen between the variations in the aromatic 
features with ionization (Fig. 1 10|) is not consistent with 
the PAH molecule model. 

There seem to be two different effects at work in HII 
regions. The first is consistent with the PAH molecule 
model and explains the behavior of the spatially resolved 
aromatic feature ratios. The second is required to ex- 
plain the good correlation between the decrease in the 
globally integrated aromatic equivalent widths and in- 
creasing ionization and the lack of a correlation between 
the aromatic feature ratios and ionization index. This 
second effect may be related to the small grain model 
for the aromatic features. If small grains can produce 
aromatic features and provide the parent population of 
the PAH mole cules in pho todissociation regions regions 
(PDRs) (|Berne et al.ll2007ft then it may be that the de- 
crease in the strength of the aromatic equivalent width is 
due to processing of the small grains. Another possible 
explanation could be that the sizes of the PDRs in re- 
gions with harder radiation fields (and, thus, higher ion- 
izations) are significantly reduced. As the aromatic fea- 
tures are known to emerge predominately from the PDR 
regions surrounding the HII regions and the underlying 
continuum emission emerges from the whole HII region, 
decreasing the area of the PDR would result in smaller 
aromatic equivalent widths. Decreasing PDR areas with 
lower metallicities may also explain the low CO emission 
seen in low metallicity starburst galaxies (jLerov et al.l 
l2007h . 

The processing of the aromatic feature carriers by mas- 
sive star formation is one line of evidence that such star 
formation actively modifies the properties of the sur- 
rounding dust grains. There is additional evidence that 
active, high-mass star formation modifies the properties 
of nearby dust grains. One observational signature of 
this effect is the lac k of the 2175 A ex t inction bump 
in sta rburst galaxies (jCalzetti et all [19941 ; I Gordon et al.l 
11997ft . Ultraviolet (UV) ex tinction curves studies i n 
the Small Magellanic Cl oud (I Gordon k. Clavtonlll998h . 
Large Magellanic Cloud dMisselt et al.lll999h. a nd Milky 
Way ([Clayton et al.ll2000t IValencic et al.ll2003l ) paint a 



consistent picture where the 2175 A bump is depressed 
and the overall UV extin ction increased in re gions near 
massive star formation (|Gordon et al.l |2003[ ). These 
changes in the UV extinction can be explained by the 
removal of the small graphite grains that produce the 
2175 A feature, along with an increase in the number of 
small silicate grains that are res ponsible for the major- 
ity of the ultraviolet ext inction ( Wcingartne r fc Draind 
1200 It IClavton et al.|[2003l ) . Given that these changes are 
in the small grain populations of the dust, it is not sur- 
prising that similar changes are seen in the mid-infrared 
where the small grains emit their absorbed energy. It is 
not clear if the changes seen in the aromatic features are 
directly correlated with the changes seen in the 2175 A 
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feature. More work in both the ultraviolet and mid- 
infrared is needed to answer this question. 

5. SUMMARY AND CONCLUSIONS 

We have studied trends in aromatic feature behavior 
with metallicity (log(0/H)+12 = 7.9-8.8) and ioniza- 
tion (II = -1.5-1.4 or [NeIII]/[NeII] = 0.03-25). Over 
this range, we have found that the equivalent width 
of the aromatic features is better correlated with the 
ionization index (II) than metallicity (log(0/H)+12) in 
both the M101 regions studied in t his paper as well 
as the starburst galaxies presented in lEngelbracht et all 
(2008). The data available for metallicities between 
log(0/H)+12 = 7.1-7.9 are not sufficient to probe these 
trends for such low metallicities. Our result implies that 
the decrease in aromatic equivalent widths seen in regions 
of massive star formation (e.g., HII regions and starburst 
galaxies) is primarily due to processing, not formation. 
This is not to say that formation does not play a role, 
but just that it is not the dominant effect. The finding of 
processing of the aromatic carriers is consistent with ex- 
isting evidence that the 2175 A extinction feature carrier 
is processed in regions of massive star formation. 

We found that the correlation of the aromatic equiva- 
lent widths with II is not a single power law as has been 
proposed in previous works. Instead the aromatic equiva- 
lent widths are constant till a threshold II of ^0 and then 
decrease. This decrease is well fit with a power law. This 



more complex beha vior is consistent with previous work 
(jBrandl et alJl2006t ) and (|Wu et al.ll2006h which probed 
a narrower range of II and seemed to find conflicting 
results. Only the large range in both II and metal- 
licity probed in this wo rk and the companion work of 
lEngelbracht et all (|2008rj allowed for the true nature of 
the aromatic feature variations to be revealed. 

The lack of a clear dependence of the ratio of different 
aromatic features to the 11.3 aromatic feature on II calls 
into question either ionization or size selective destruc- 
tion of PAH molecules as the mechanism for the observed 
reduction in the aromatic equivalent width versus ion- 
ization. This mig ht indicate a more solid material origin 
for the aromatics (|Berne et al.ll2007T ) which dominate the 
global measurements of HII regions. Or it might imply 
that the area of the PDR regions decreases with respect 
to the area of the small grain emission as the radiation 
field hardness (and thus ionization) increases. 
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TABLE 4 

Atomic Emission Line Strengths 3 - 



name 


[Aril] 


[Arlll] 


[SIV] 


[Nell] 


[NcIII] 


[SIII] 18 


[SIII] 33 


[Sill] 




7.0 fim 


9.0 fim 


10.5 fim 


12.8 fim. 


15.5 fim 


18.7 fim 


33.5 /im 


34.8 fj,m 


Nucleus 


3.30 ± 0.51 






7.57 ±0.51 


1.20 ±0.25 


4.50 ±0.51 


18.78 ±0.89 


20.26 ± 0.79 


Hodge 602 


1.00 ±0.23 




0.34 ±0.12 


2.30 ±0.16 




1.54 ±0.13 


3.74 ±0.12 


3.48 ±0.14 


Searle 5 


1.04 ±0.17 


0.09 ± 0.09 


0.18 ±0.05 


4.13 ±0.16 


0.16 ±0.11 


4.01 ±0.18 


20.36 ±0.28 


9.14 ±0.18 


NGC 5461 


1.78 ± 5.33 


10.74 ±0.81 


28.39 ±0.81 


27.32 ± 1.54 


56.09 ±5.82 


60.56 ±3.94 


259.21 ± 5.53 


80.58 ± 3.67 


NGC 5462 


0.86 ± 0.53 


2.25 ± 0.50 


2.75 ± 1.21 


4.34 ± 7.83 


6.50 ±1.95 


7.46 ± 1.49 


43.90 ±1.02 


15.01 ± 0.70 


NGC 5455 


0.34 ± 1.61 


1.25 ±0.09 


3.08 ± 0.09 


3.85 ±0.15 


11.19 ± 3.19 


7.93 ± 1.68 


46.14 ±0.74 


18.63 ±0.59 


NGC 5447 


0.25 ± 0.11 


0.92 ± 0.11 


3.33 ±0.10 


1.61 ±0.20 


4.27 ±0.56 


5.41 ±0.64 


27.51 ± 1.31 


10.60 ± 1.07 


NGC 5471 




0.52 ± 0.15 


12.92 ±0.73 


0.96 ±0.29 


20.63 ± 2.33 


8.47 ± 1.38 


32.60 ±0.75 


11.27 ±0.57 



3 Units are TO -17 W rrT 2 . 



TABLE 5 
Aromatic Emission Line Strengths 3 



name 5.7 (jm 6.2 fim 7.7 (im 8.3 fim 8.6 (im 10.7 (im 11.3 fim 12.0 fim 12.7 (im 14.0 (im 17.0 (im 



Nucleus 7.86 ± 1.87 77.06 ± 1.78 271.43 ±7.87 23.36 ± 2.00 56.62 ±2.14 1.64 ± 0.47 61.58 ± 2.36 15.98 ± 1.07 36.27 ± 2.27 1.04 ±0.57 55.12 ±3.72 

Hodge 602 2.53 ± 1.25 9.61 ±1.16 37.42 ± 3.83 1.76 ± 0.90 16.03 ± 0.95 1.37 ±0.24 12.08 ±1.03 1.63 ± 0.38 4.25 ± 0.74 ■■• 4.97 ± 0.42 

Searle 5 •■■ 12.54 ± 0.44 41.64 ± 2.20 1.03 ± 0.49 9.77 ± 0.45 0.02 ±0.13 9.39 ± 0.45 1.88 ±0.22 4.73 ± 0.50 0.59 ±0.14 1.99 ±0.30 £ 

NGC 5461 6.64 ±1.25 66.21 ± 2.58 205.47 ± 13.74 12.74 ± 3.50 40.30 ± 3.08 2.48 ± 1.01 50.56 ±4.50 14.73 ± 2.41 30.49 ± 5.56 ■■• 3.21 ±4.16 § 

NGC 5462 1.46 ±1.14 20.21 ± 1.54 64.84 ± 7.92 4.42 ± 1.86 13.38 ± 1.89 1.34 ± 0.64 16.13 ±2.42 1.88 ±1.30 13.54 ± 3.05 2.95 ± 0.71 0.72 ± 1.18 o 

NGC 5455 0.78 ±0.64 14.42 ± 0.48 38.75 ± 2.86 2.84 ± 0.49 10.77 ±0.42 0.39 ±0.16 10.87 ±0.51 2.34 ± 0.32 5.93 ± 0.65 0.63 ±0.21 ■■• P 

NGC 5447 0.50 ±0.44 1.30 ±0.52 8.41 ± 1.58 1.35 ± 0.74 2.12 ±0.68 ■•■ 2.93 ± 0.75 0.54 ±0.37 1.90 ±0.78 ■■• 2.67 ± 1.24 g. 

NGC 5471 0.24 ±0.69 •■■ 8.82 ± 0.76 ■■• 2.15 ±0.46 1.27 ±0.29 2.85 ± 0.92 •■■ 0.84 ± 0.38 0.74 ± 0.49 0.58 ± 1.73 m 



a Units arc 1CT 17 W m~ 2 . 



TABLE 6 

Aromatic Emission Line Equivalent Widths 3 



name 5.7 (im 6.2 fim 7.7 fim 8.3 fim 8.6 fim 10.7 fim 11.3 fim 12.0 fim 12.7 fim 14.0 fim 17.0 fim 



Nucleus 0.16 ±0.04 1.44 ±0.03 4.91 ±0.14 0.46 ± 0.04 1.17 ±0.04 0.05 ± 0.01 1.90 ±0.07 0.51 ± 0.03 1.18 ±0.07 0.03 ± 0.02 1.46 ±0.10 

Hodge 602 0.29 ±0.14 1.21 ±0.15 6.09 ± 0.62 0.35 ±0.18 3.80 ± 0.23 0.54 ± 0.09 3.92 ± 0.34 0.41 ±0.10 0.89 ±0.16 ••■ 0.59 ± 0.05 

Searle 5 ••■ 1.47 ±0.05 4.85 ± 0.26 0.13 ±0.06 1.36 ± 0.06 0.00 ± 0.03 1.81 ± 0.09 0.33 ± 0.04 0.78 ±0.08 0.08 ± 0.02 0.17 ±0.03 

NGC 5461 0.13 ±0.02 0.99 ± 0.04 2.07 ±0.14 0.12 ±0.03 0.36 ± 0.03 0.02 ± 0.01 0.34 ± 0.03 0.09 ± 0.01 0.16 ±0.03 ••■ 0.01 ± 0.01 

NGC 5462 0.11 ±0.08 1.36 ±0.10 4.35 ± 0.53 0.32 ±0.13 0.99 ±0.14 0.12 ± 0.06 1.37 ±0.21 0.15 ±0.10 0.95 ±0.21 0.15 ±0.04 0.02 ± 0.03 

NGC 5455 0.07 ±0.06 1.12 ±0.04 2.83 ± 0.21 0.21 ± 0.04 0.82 ± 0.03 0.03 ± 0.01 0.75 ± 0.04 0.14 ±0.02 0.31 ± 0.03 0.02 ± 0.01 

NGC 5447 0.14 ±0.12 0.29 ±0.12 1.31 ±0.25 0.19 ±0.11 0.29 ± 0.09 •■• 0.30 ± 0.08 0.05 ± 0.03 0.16 ±0.06 ■•■ 0.10 ±0.04 

NGC 5471 0.08 ±0.22 ••■ 1.60 ±0.14 ■•• 0.32 ± 0.07 0.11 ±0.03 0.21 ± 0.07 •■■ 0.04 ± 0.02 0.03 ± 0.02 0.01 ± 0.03 



Units /xm. 



